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I describe CLEO-c purely leptonic decay results leading to / 
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D+ = (222.6 ± 16.7±1;|) MeV, 
/ c+ = (280.1 ± 11.6 ± 6.0) MeV, and f D +/f D + = 1.26 ± 0.11 ± 0.03. Form-factor measurements 
in Cabibbo favored and suppressed pseudoscalar decays are presented. Some comparisons are made 
with theoretical predictions. 



I. INTRODUCTION 



Threshold production of D°D and D + D mesons at 
3770 MeV, and D+D*- + D*+ D~ mesons at 4170 MeV 
in e + e~ annihilations have allowed CLEO-c to make pre- 
cision measurements using purely leptonic and semilep- 
tonic charm meson decays. 



II. PURELY LEPTONIC DECAYS 

To extract precise information from B mixing mea- 
surements the ratio of "leptonic decay constants," fi for 
Bd and B s mesons must be well known [l|. Indeed, 
the recent measurement of B2 mixing by CDF @ has 
pointed out the urgent need for precise numbers. The fi 
have been calculated theoretically. The most promising 
of these calculations are based on lattice-gauge theory 
that include the light quark loops [|[. In order to ensure 
that these theories can adequately predict /s s //s d it is 
critical to check the analogous ratio from charm decays 
f D + / fn+ ■ Here I present the most precise measurements 
to date of f D +, f D + 0,0 and f D +/f D +. 




FIG. 1: The decay diagram for 



In the Standard Model (SM) the LXK meson decays 
purely leptonically as shown in Fig. [1] The decay width 
is given by @ 



r(£>+ = ^f n+ m 2 e+ M 



G> 2 
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(1) 



x U-m 2 i+ /M 2 D +) \V cq \ 2 , 



where m«+ and M n + are the I 



and D+ s) 



masses, \V cq \ 



is the CKM element appropriate to either D + (V c d) or 
Df (V cs ) decay and Gp is the Fermi constant. 

New physics can affect the expected widths; any undis- 
covered charged bosons would interfere with the SM W + . 
These effects may be difficult to ascertain, since they 
would simply change the value of the fi's. The ratio 
f D + 1 fry+ , however, is much better predicted in the SM 
than the values individually. Akeroyd predicts that the 
presence of a charged Higgs boson would suppress this 
ratio significantly In addition, the ratio of decay 
rates to different leptons are fixed only by well-known 
masses in Eq. [T] For example, the SM prediction for 
r(Z>+ -> r+i/)/r(L>+ -> n+v) is 9.72. In general, any 
deviation from a predicted ratio would be a manifesta- 
tion of physics beyond the SM, and would be a clear 
violation of lepton universality [8( . 

CLEO previously measured f D + using 4.8 fb -1 of con- 
tinuum annihilation data at or just below the T(4S) [||. 
This analysis introduced a number of new ideas: (i) The 
7 and fi + from D* + — > r )D + \ Df — * fi + i> were detected 
directly, and the v 4- vector was inferred from missing en- 
ergy and momentum measurement in half of the event, 
where the event half was determined using the normal 
to the thrust axis, (ii) The v 4-vector was corrected to 
get the right D + mass and AM = M(^[± + v) — M(fi + is) 
was examined (see Fig. EJa)). (ii) The background was 
measured using the same technique with e + identified 
instead of , relying on the large suppression of the 
e + rate compared with the fi + rate, (iv) The reaction 
D*° -> 7D ; D° -> K--K+, where the tt+ is first found 
and then ignored was used to evaluate efficiencies. The 
published result was 



r(D+ - fj+v) 



t(d; 



= 0.173 ± 0.023 ±0.035 



(2) 



BaBar recently performed an improved analysis based 
on these techniques [l0|. They used 230 fb -1 of con- 
tinuum data. To reduce the background and systematic 
errors they fully reconstruct a D°, D + or D* meson in 
the event with the 7 and fi + candidate. Their data are 
shown in Fig. (Hb) . They find 
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r(D+-/i+i/) 
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= 0.143 ±0.018 ±0.006 



(3) 



2 




300 






r (b) rt 


250 






200 






150 




I 


100 


1 

. JJl / 




50 




4UlJB iff 








0.05 0.1 



0.15 0.2 0.25 
AM [GeV/c 2 ] 



FIG. 2: The AM distributions for /i + candidates after the e + 
subtraction from (a) CLEO and (b) BaBar. The solid curves 
are fits to signal plus background. 



Both of these results, however, need to assume a value 
for B(D+ — > 07r + ) [ll], in order to extract the decay 
constant. Because of interferences among the final state 
K + K~ir + particles, the rate for </>7r + depends on exper- 
imental cuts fl2l |. and thus has an inherent, sizable, sys- 
tematic error. (Other experiments also normalize with 
respect to this or other less well known modes.) 

CLEO-c eliminates this uncertainty by making abso- 
lute measurements. We tag a D~ decay and search for 
three separate decay modes of the D+: (1) and t + v, 



where (2) r + — > tt + V or (3) r + — > e + z47 [l3j. For the first 
two analyses we require the detection of the 7 from the 
D* — > "fD s decay, irrespective if the D* a is the parent 
of the tag or the leptonic decay. In either case, for real 
D*D S events, the missing mass squared recoiling against 
the photon and the Dj tag should peak at M D + and is 
given by 

MM* 2 = (£ CM -E D - E y f - ( p> C M - ~Pd - ~P-,f , 

where Eqm (Pcm) is the center of mass energy (momen- 
tum), Ed (pd) and E 1 (jp?j) are the energy of the fully 
reconstructed D~ tag, and the additional photon. In 
performing this calculation we use a kinematic fit that 
constrains the decay products of the D~ to M D + and 
conserves overall momentum and energy. 

The MM* 2 from the D~ tag sample data is shown in 
Fig. El There are 11880±399±511 signal events in the 
interval 3.978 >MM* 2 > 3.776 GeV 2 . 
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FIG. 3: The MM* 2 distribution from events with a photon in 
addition to the D~ tag. The curve is a fit to the Crystal Ball 
function and a 5th order Chebychev background function. 



Candidate Df — > n + v events are searched for by se- 
lecting events with only a single extra track with opposite 
sign of charge to the tag; we also require that there not 
be an extra neutral energy cluster in excess of 300 MeV. 
Since here we are searching for events where there is a 
single missing neutrino, the missing mass squared, MM 2 , 
evaluated by taking into account the seen fi + , Dj, and 
the 7 should peak at zero, and is given by 



MM 2 = (E CM - E D - E 1 - E^ 
- (~Pcm - ~Vn-~$i- 



(4) 



P, 



where E^ (p/j,) is the energy (momentum) of the candi- 
date muon track. 

We also make use of a set of kinematical constraints 
and fit the MM 2 for each 7 candidate to two hypotheses 
one of which is that the D~ tag is the daughter of a D* ~ 
and the other that the D* + decays into jDf, with the 
Df subsequently decaying into [i + v. 

The kinematical constraints are the total momentum 
and energy, the energy of the either the D* or the D s , 
the appropriate D* — D s mass difference and the invariant 
mass of the D s tag decay products. This gives us a total 
of 7 constraints. The missing neutrino four-vector needs 
to be determined, so we are left with a three-constraint 
fit. We perform a standard iterative fit minimizing % 2 . 
As we do not want to be subject to systematic uncertain- 
ties that depend on understanding the absolute scale of 
the errors, we do not make a x 2 cu t, but simply choose 
the photon and the decay sequence in each event with 
the minimum \ 2 ■ 

We consider three mutually exclusive cases: (i) the 
track deposits < 300 MeV in the calorimeter, charac- 
teristic of a non-interacting 7r + or a (ii) the track 
deposits > 300 MeV in the calorimeter, characteristic of 
an interacting 7r + ; (iii) the track satisfies our e + selection 
criteria. The MM 2 distributions are shown in Fig. 2] The 
separation between fi + and tt + is not unique. Case (i) 
contains 99% of the /i + but also 60% of the tt + , while 
case (ii) includes 1% of the /i + and 40% of the ir + [f|. 
There is a clear peak in Fig.HJi), due to Df — > Fur- 
thermore, the events in the region between peak and 
0.20 GeV 2 are dominantly due to the t + v, t + — > n + v 
decay. The best result comes from summing case (i) and 
case (ii) below MM 2 of 0.20 GeV 2 ; higher values of MM 2 
admit background from r]ir + and K°ir + final states. The 
branching fractions are summarized in Table |TJ The ab- 
sence of any detected e + opposite to our tags allows us 
to set the upper limit listed in Table [U 



CLEO-c also uses Dt 



,+ 



vv. Electrons of 



opposite sign to the tag are detected in events without 
any additional charged tracks, and determining the un- 
matched energy in the crystal calorimeter (EgQ ra ). This 
energy distribution is shown in Fig. [5] Requiring E 1 



cc 



< 



400 MeV, enhances the signal. The branching ratio re- 
sulting from this analysis is also listed in Table |TJ 

CLEO-c's published result for f D + Q uses the "double- 
tag" method at 3770 GeV, where D + D~ final states are 
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FIG. 4: The MM 2 distributions from data using D~ tags 
and one additional opposite-sign charged track and no extra 
energetic showers (see text). 
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FIG. 5: The extra calorimeter energy from data (points), com- 
pared with the Monte Carlo simulated estimates of semilep- 
tonic decays in general (dotted), the K°e + v mode specifically 
(shaded), as a sub-set of the semileptonics, and the expecta- 
tion from signal (dashed). The peak near 150 MeV is due to 
the 7 from —* ^yDs decay. (The sum is also shown (line).) 
The arrow indicates the selected signal region below 0.4 GeV. 



TABLE I: Measured D + Branching Fractions 



Final State 



fj, + v 



g (%) 



(r+ -> tt+v) 
(t + — > e + uv) 



tv (average) 



0.657 ±0.090 ±0.028 
0.664 ± 0.076 ± 0.028 
7.1 ± 1.4 ±0.3 
6.29 ±0.78 ±0.52 
6.5 ±0.8 
< 3.1 x 10~ 4 (90% cl) 



t From summing the fi + u and t + v contributions for MM 2 < 
0.20 GeV 2 . 



produced without any extra particles. Here one D~ is 
fully reconstructed and then there are enough kinematic 
constraints to search for D + — > fi + v by constructing the 
missing mass-squared (MM 2 ) opposite the D~ and the 
muon. Fifty signal events are found of which 2.8 are 
estimated background, resulting in: 



= (4.40 ± 0.66 



+0 of 2 ) x 10- 



(5) 



The decay constant f D + is obtained from Eq. (TT|) us- 
ing 1.040±0.007 ps as the D+ lifetime, and \V cd \ = 
0.2238±0.0029, giving 



h 



(222.6 ± 16. 7±§;f) MeV 



(6) 



CLEO-c also sets limits on B(D + — > e + v e ) < 2.4 x 
10^ 5 , [1] and B(D + — > t + v) branching ratio to < 2.1 x 
10~ 3 at 90% CL. [ll|- These limits are consistent with 
SM expectations. 

For Dj decays, we first test lepton universality in 



consistent with the predicted value of 9.72. Combining 
our branching ratios determinations and using r D +=0.49 
ps and |V CS |=0.9737, we find 

f Da = (280.1 ± 11.6 ± 6.0) MeV, and (8) 
f D + //d+ = 1-26 ± 0.11 ± 0.03. 

These preliminary results are consistent with most re- 
cent theoretical models. As examples, unquenched lattice 
[H predicts 1.24 ± 0.01 ± 0.07, while one quenched lat- 
tice calculation [16| gives 1.13 ± 0.03 ± 0.05, with other 
groups having similar predictions [201 ] . 



III. SEMILEPTONIC DECAYS 

One of the best ways to measure magnitudes of CKM 
elements is to use semileptonic decays since they are far 
simpler to understand than hadronic decays and the de- 
cay width is ~ iV'cql 2 . On the other hand, measurements 
using other techniques have obtained useful values for 
V cs and V c d [13] > and thus semileptonic D decay mea- 
surements are a good laboratory for testing theories of 
QCD. For a D meson decaying into a single hadron (h), 
the decay rate can be written exactly in terms of the 
four-momentum transfer defined as: 



q = (p£ - Ph) = m D+ m h- 2E h m D 



(9) 



For decays to pseudoscalar mesons and "virtually 
massless" leptons, the decay width is given by: 



T+V) 



r{D+ -» fi+v) 



= 9.9 ± 1.9 



(7) 



dT(D -> Pe+i 



\V cq \ 2 Glp 3 P 

247T 3 



|/+(9 2 )| , 



(10) 



4 



where pp is the three-momentum of P in the D rest 
frame, and f+(q 2 ) is a "form- factor," whose normaliza- 
tion must be calculated theoretically, although its shape 
can be measured. 

The shape measurements can distinguish between 
form-factor parameterizations. In general, 



Uii 2 ) = 



U(o) 



1 



(l-Op)(l 



poll 



) 



(Md+Mp) 2 



/2 Imf(q' 2 ) 



which incorporates the possibility of a virtual of a nearby 
pole (first term) with fractional strength a p . The integral 
term can be expressed in terms of an infinite series [18j |. 
Typically it takes only a few terms to describe the data. 
An analytical parametrization 



U(q 2 



/+(o) 



(1 - q 2 /m 2 pole )(l - aq 2 /m 2 polo ) 



(11) 



has become popular (l9( , though it has been criticized as 
being overly constraining [l8| . Fits are typically done for 
/+(0) and either m po i c or a. Naively, setting a to zero 
gives the simple pole model where the pole mass corre- 
sponds to the first vector resonance in the D-P system, 
D* s for D -> Ke~u and D* for D — » ixe~v 

CLEO-c uses two methods to analyze pseudoscalar de- 
cays. In the first method tags are fully reconstructed and 
events with a missing v are inferred using the variable 
U = -Emiss — |imiss|) similar to MM 2 , where "miss" here 
refers to the missing energy or momentum (see Fig. [6]). 
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FIG. 6: U distributions using D tags in conjunction with an 
identified electron of opposite flavor plus a single hadron. The 
peak centered at zero is signal. The dashed curves indicate 
various backgrounds, while he solid curve shows the fit to 
signal plus background. 
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FIG. 7: Mbc distributions for events containing an identified 
electron plus a single hadron candidate. The shaded regions 
indicate various backgrounds. 
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FIG. 8: CLEO-c form-factor shapes using the tagged sample. 
The lower curves are fits to the modified pole model, while 
the upper curves are fits to unquenched lattice QCD [2ll |. 



significant disagreement. Properties of these decays are 
listed in Table GEO 

Measurements of the vector decays D — > K*e + v and 
pe + u can be used to determine \V u b\ along with measure- 
ments of B -» pl~v and B -> K*£ + £- 0. CLEO-c has 
examined D vector semileptonic decays. Non-parametric 
form-factors in the Cabibbo favor D° — ► K*°e + v de- 
cays have been measured by CLEO-c (24|, following a 
method developed by FOCUS [25j . Cabibbo suppressed 
form-factors have been measured in D — > pe + v de- 
cays. The U distribution for pe + v decays is shown in 
Fig. [5J Preliminary branching fractions are listed in Ta- 
ble IIIII along with observations (or limits) from other 
rare semileptonic decays. Selected candidates are used 
to measure the ratios of pole dominated form-factor ra- 
tios as R v = 1.40±0.25±0.03 and R 2 = 0.57±0.18±0.06, 



The second method consists of also using missing en- 
ergy and momentum, skipping the step of reconstructing 
the tag, but using all of the measured charged tracks and 
photons. Then the D mass is reconstructed. The beam- 
constrained mass (Mbc) distributions are shown in Fig. [7] 

Both cases have excellent signal to background in these 
modes. The ^-reconstruction has better statistical albeit 
poorer systematic errors. Eventually combined results 
will be quoted; they should not be averaged as there are 
a substantial number of events in common. 

Form-factor shapes using the tagged sample are are 
shown in Fig. [51 The unquenched lattice QCD model 
PH is systematically higher than our data, but not in 



TABLE II: Properties of D° — ► P e + v decays (preliminary) 
[221 ] . To determine m po i e , a in Eq. [TTJ is set to zero. 



Quantity 


K~e+v 


7T e V 


Source 




B{%) 


3.58(5)(5) 


0.309(12)(6) 


CLEO-c 


Tag 


B{%) 


3.56(3)(11) 


0.301(11)(10) 


CLEO-c 


NoTag 


B{%) 


3.58(18) 


0.360(60) 


PDG04 




l/+(0)| 


0.761(10)(7) 


0.660(28)(11) 


CLEO-c 
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1.96(3)(1) 


1.95(4)(2) 


CLEO-c 


Tag 


ra po i c (GeV) 


1.97(2)(1) 


1.89(3)(1) 


CLEO-c 


NoTag 


a 


0.22(5)(2) 


0.17(10)(5) 


CLEO-c 


Tag 


a 


0.21(4) (3) 


0.32(7) (3) 


CLEO-c 


Tag 



5 




U (GeV) U (GeV) 



FIG. 9: U distributions for (a) D + and (b) D° decays into 
pe + v candidates. The dashed curve shows the signal, the 
dotted curves show various backgrounds and the solid curve 
the sum. 



TABLE III: Rare semileptonic decay branching fractions 



Decay 


B x 10~ 4 


D" -> p-e+v 


15.6 ±1.6 ±0.9 


D+ -» p°e+u 


23.2 ±2.0 ±1.2 


D+ -> Loe + u 


14.9 ±2.7 ±0.5 


D+ -> (pe+v 


< 2 at 90% CL 


D + — > ?7e + !/ 


12.9 ±1.9 ±0.7 


D+ -> r/e+f 


< 3 at 90% CL 


D" -> K-TT + -n-e + v 


2.9lio ±°- 5 



using both charge and neutral modes [22| . 

Other results on semileptonic decays from CLEO-c in- 
clude measurement of the inclusive D° and D + semilep- 
tonic branching fractions of (6.46 ± 0.17 ± 0.13)%, and 
(16.13 ± 0.20 ± 0.33)%, respectively, leading to a mea- 



surement of the partial width ratio of T(D + )/T(D°) = 
(0.985 ±0.028 ±0.015), consistent with isospin symmetry 

IV. CONCLUSIONS 

CLEO-c measurements of leptonic and semileptonic 
decays have already reached precisions that provide very 
useful benchmarks for testing of QCD theories. From 
leptonic decays we have 

f D + = (222.6 ± 16.7±1;|) MeV, (12) 
f D + = (280.1 ±11.6 ±6.0) MeV, 
f D +/fo+ = 1.26 ± 0.11 ± 0.03 . 

These results are consistent with most theoretical calcu- 
lations including those of unquenched lattice QCD (20j . 

CLEO-c is also breaking new ground in the study of 
semileptonic decays. Form-factors in Cabibbo suppressed 
decays are reaching an unprecedented level of accuracy 
and are also confronting theory. 
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